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Abstract 

We treat the problem of the approach to thermal equilibrium by only resorting 
to quantum dynamics of an isolated macroscopic system. Inspired by the two 
important works in 2009 and in 1929, we have noted that a condition we call 
"thermodynamic normality" for a macroscopic observable guarantees the approach 
to equilibrium (in the sense that a measurement of the observable at time t almost 
certainly yields a result close to the corresponding microcanonical average for a 
sufficiently long and typical t). A crucial point is that we make no assumptions 
on the initial state of the system, except that its energy is distributed close to a 
certain macroscopic value. 

We also present three (rather artificial) models in which the thermodynamic 
normality can be established, thus providing concrete examples in which the ap- 
proach to equilibrium is rigorously justified. 

Note that this kind of results which hold for any initial state are never possible 
in classical systems. We are thus dealing with a mechanism which is peculiar to 
quantum systems. 

The present note is written in a self-contained (and hopefully readable) man- 
ner. It only requires basic knowledge in quantum physics and equilibrium statis- 
tical mechanics. 
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1 Background and main results 

1.1 Setup, background, and motivation 
Setup 

Consider a finite but macroscopic quantum mechanical system which is completely iso- 
lated from the outside world. The system is fully described by its Hamiltonian H. For 
a = 1, 2, . . ., we denote by and ipa the eigenvalue and the corresponding normalized 
eigenstate of H. We assume that the eigenvalues are non-degenerate, i.e., 7^ Ea' 
whenever a 7^ a' . 

We want to discuss the thermal equilibrium state of the system when it has a macro- 
scopic energy E. The standard (and empirically justified) procedure is to invoke the 
"principle of equal probability", and declare that all the microscopic states within the 
energy rang^ [E^ E + AE] contribute equally to the equilibrium state. Here AE is a 
small, but still macroscopic energy. 

To be precise let A be any (macroscopic) observable of the system. Its microcanonical 
average at energy E is defined bjfl 

{Ea(i[E,E+AE]) {Ea&[E,E+AE]) 

where we denoted by He,ae the Hilbert space spanned by all the energy eigenstates ipa 
such that Ea G [E,E + AE], and Tr^^ ^bI' ' ' ] denotes the trace over this space. 

^ The choice of this range is rather arbitrary. The only requirements are that (i) macroscopic 
properties of the system remain almost the same through the range, and (ii) there are huge number 
of energy eigenstates within the range. It is misleading to imagine that AE is determined by physical 
processes such as the measurement of the energy. 

When A is defined in terms of B, we write A := B ot B —: A. 
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Our goal is to justify the use of the microcanonical average (11. ip by resorting only 
to quantum mechanics. 

Initial state and time evolution 

Suppose that at initial time t = 0, the system was in a quantum mechanical pure state 
f{0). We assume that the initial pure state has energy distributed near E. For simplicity 
we let y^{0) be an arbitrary normalized state in the above defined space 'He,ae- Of course 
'He,ae contains many states which are radically far from equilibrium. 
Now write this initial state as 

¥^(0) = J]c«V^a (1.2) 

a 

where we have Ca = Q Ea ^ [E^E + AE], and Ylia 1*^"^ ~ '^^^ time evolution of 
this state is given by 

<^W = $^c„e-^^"V„ (1.3) 

a 

for any t > 0. Our expectation is that for sufficiently long and typical time t, macroscopic 
properties of the pure state (pit) can be described by the microcanonical average (II. ip 
in the sense that a measurement of A almost certainly yields a result close to {A)^c- 
This expectation is based on the standpoint called "individualist" point of view in [Ij. 
See [1] for a list of related work^. 

Basic picture 

Let us briefly discuss the idea behind the above expectation that ip{t) should be described 
by the microcanonical average. See [2j (and [3j, if the reader can read Japanese) for 
more details. 

We first argue that the universal applicability of equilibrium statistical mechanics 
strongly suggests that (i) overwhelming majority of sates in the space 'He,ae look almost 
similar from macroscopic points of view, and (ii) "equilibrium properties" are nothing 
but the common properties that are shared by these majority of states. In other words, 
most of the states in 1-Le,ae can be regarded (from macroscopic points of view) as rep- 
resentatives of the equilibrium state. The remaining exceptional minorities correspond 
to varieties of nonequilibrium states. See Figure [1] It is then evident that by averaging 
over all the states in 'He,ae as in (II. ip . one can extract the properties of the majority 
since exceptional contributions are so minor and averaged out. 

Once accepting this picture, we can naturally see why (pit), in the long run, should 
approach equilibrium. Suppose that the initial state ipiQi) describes a physical situation 

^ It is sometimes argued that one should consider a system coupled to the outside world because no 
systems are perfectly isolated in reality. Let us emphasize, however, that it is perfectly reasonable to 
consider an isolated system as an idealization. After convincing ourselves that an isolated macroscopic 
system has a tendency to reach its equilibrium, we can start thinking about why and how the coupling 
to the outside world does not destroy this tendency. 
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Figure 1: The big elliptic region represents the Hilbert space T-Le,ae, the space 
of states that have a macroscopic energy E. Overwhelming majority of the states 
look almost the same from macroscopic points of view, and they exhibit properties 
of the corresponding equilibrium states. Exceptional states which belong to very 
small shaded regions have macroscopic properties different from equilibrium. Note 
that this is a very crude picture of the Hilbert space with extremely high dimension. 
Also note that the exceptional regions must be much smaller than depicted. The 
figure is taken from [3]. 




Figure 2: The basic mechanism of the approach to equilibrium based on the 
picture as in Figure [TJ Left: If the state was initially in one of the exceptional 
regions, it eventually moves out of the region and evolves into the overwhelming 
majority of typical states. Right: If the state was initially typical, it rarely wanders 
into the exceptional regions, and keeps being typical. The figure is taken from [3]. 
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which is very far from equihbrium, e.g., two bodies at different temperatures in contact 
with each otheil§. The initial state V'(O) certainly is an "exceptional" state which does 
not belong to the overwhelming majority. Obviously the coefficients in fll.2p must 
be chosen with an extreme care and precision in order to realized such an exceptional 
state. As t grows, each coefficient gets individual phase factor as in (11. 3p . and the very 
delicate balance of the coefficients realized at t = will soon be lost. It is quite likely 
that, after a sufficiently long time, the state ip(t) is no longer exceptional, and belongs 
to the overwhelming majority of equilibrium states. See Figure [2J 

Older wor 

If one wants to justify the above expectation that ip{t) becomes an equilibrium state, a 
natural (and well-known) starting point is to consider a long-time average. 

Let A be a macroscopic observable. Its quantum mechanical expectation value in 
the state ip(t) is 

((p(t),i(^(t)> = J]c^c„,e^(^"-^'^')*(V^„,i^„,), (1.4) 

a,a' 

where we used (II. 3p . Since we assumed that the energy eigenvalues are nondegenerate, 
all the terms with a ^ a' exhibit oscillations. By averaging over infinite amount of time, 
the oscillations simply average out, and one finds 

limi/ dt{if{t),Aip{t)) = y2\c^mi;^,A^^), (1.5) 

" a 

which is remarkably similar to the microcanonical average (II. ip . 

Then the challenge was to show that the quantum mechanical expectation value 
(without the time average) is close to the microcanonical average. More precisely, to 
show that 

{ip{t),Aip{t)) ^ {A)mc for most t (1.6) 

provided that we consider long enough time interval. This, if done, may be regarded as 
a justification of the microcanonical averag^. 

In the temporal fluctuation of (^ip{t), Atpit)'^ was evaluated using the standard 
Chebyshev inequality argumentl^. It was shown that the desired (11.61) is justified under 
the assumption that each coefficient Cq is smal0 In other words, the initial state V'(O) 



^ We treat (an artificial version of) this problem in section 1231 

^ Note that this part of the note is a rather personal account which places the present observation 
in a broader context. Please refer to pLj and references therein for a more balanced view on what have 
been done. 

^ To be more precise, one further has to assume that quantum mechanical fluctuation of A in Lpit) 
is small. 

^ It was also assumed that ~ E^' — Ea" — Ea'" ^ implies a = a" , a' = a'" . This is called the 
non-resonance condition. 

[4] discusses a setup leading to the canonical ensemble, but the microcanonical setting is treated 
in exactly the same (indeed easier) manner. See also [5] where we treated a microcanonical problem. 
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should have energy distributed around E, but not too sharply. A different criterion of 
the same nature was recently obtained in [B]. See section [Q] for a new result of the 
same philosophy. 

It is a very difficult question whether we should be satisfied with this type of "deriva- 
tion of statistical mechanics from quantum mechanics." The most delicate point is the 
restriction on the allowed initial state v^(0). It seems to be true that excluded initial 
states (which have sharp peaks at some energies) are "rare" , and the practical "chance" 
of choosing such initial states in actual world can be safely neglected. But in what sense 
are they "rare" or do they have negligible "chances" ? We are dealing with a preparation 
of a pure state, where we do not see any natural concept of probabilities. 

A new (but in fact rather old) thought 

In [7], Goldstein, Lebowitz, Mastrodonato, Tumulka, and Zanghi proposed an abstract 
formulation of equilibrium states in a macroscopic system. They were able to prove 
that, in a typical system, any initial state ^^(O) G He,ae will eventually evolve into 
equilibrium. In other words the approach to equilibrium is guaranteed by only assuming 
that the initial state has energy distributed around E. Note that this assumption on 
the range of energy is almost mandatory since we never observe a linear combination of 
multiple states with macroscopically different energieJ^. 

Such a result with basically no restrictions on the initial state has a great advantage, 
at least from a conceptual point of view. With such a statement, we no longer have 
to worry about the "chance" of choosing an exceptional initial state. It should be 
stressed at this moment that, in a classical system, it is impossible to have a similar 
statement which holds for any initial states. In any classical dynamical system which 
shows thermalization, there exists a set of "exceptional" initial states which do not 
relax to equilibrium. The set of exceptional states usually have a vanishing (Lebesgue) 
measure, and one never chooses such a state provided that the initial state is sampled 
according to the microcanonical ensemblj^. 

Interestingly (and surprisingly), back in 1929, a result of the same nature as that 
of Goldstein, Lebowitz, Mastrodonato, Tumulka, and Zanghi was already derived by 
von Neumann [S]. He treated the problem of justifying statistical mechanics by using 
quantum mechanics, and stated a theorem which holds for any initial states. We highly 
recommend the readers to go through [1], which is a very well written commentary by 
Goldstein, Lebowitz, Tumulka, and Zanghi on this important (and long forgotten) work 
by von Neumann. 

We will discuss the result of Goldstein, Lebowitz, Mastrodonato, Tumulka, and 

It is possible that such a state (Schrodinger's cat) is inhibited (or destroyed) by interaction of the 
system with the outside world. But this is a totally different topic (for the moment). 

Although this is usually regarded as a satisfactory justification of the equilibrium statistical me- 
chanics, it should be noted that here one is invoking a kind of circular reasoning. There seems to be no 
a priori reason (based solely on classical mechanics) to use the microcanonical (or the Lebesgue) mea- 
sure when one samples the initial state. It seems possible that the uncertainty principle from quantum 
mechanics provides a reason. 
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Zanghi [7] agin in section 11.41 



1.2 Definition of "thermodynamic normality" and the main 
theorem 

We are now ready to describe our simple observation inspired by the results of Goldstein, 
Lebowitz, Mastrodonato, Tumulka, and Zanghi [7] and also of von Neumann [8J. Like 
these authors, we insist on deriving a result that is valid for any choice of initial state 
ip{0) whose energy is distributed around E. We use the same notation as in section [LTl 
The key of the observation is the following definition of thermodynamic normality of 
a macroscopic observablj^. 

Definition: We say that a macroscopic observable A is thermodynamically normal with 
respect to H in the interval [E, E + AE] if 

(V'„,(i-(i)„,c)Va> <^'C (1-7) 

holds for any a such that E^ G [E,E + AE], where C > is a small constant. Here 
constant A (which may depend on E) is a (rather arbitrarily chosen) typical magnitude 
of i. 

Note that by averaging the above fll.7p over all a such that E^ G [E,E + AE], one 
gets 

{{A-{AU^r)^^^<A\ (1.8) 

which simply says that the fluctuation of A is small in the microcanonical ensemble. 
This property is believed to hold quite generally, and may be proved rigorously under 
suitable concrete settings (although the proof may not be easy). 

By the thermodynamic normality (11.71) . on the other hand, we are requiring that, as 
far as the observable A is concerned, each eigenstate ipa in the energy range [E, E + AE] 
behaves almost as the equilibrium state. This is indeed a very strong requirement, which 
may or may not be true for general macroscopic systems. In section [2], we discuss (rather 
artificial) examples where the thermodynamic normality can be justified. See section [3] 
for further discussions. 

The following main theorem is indeed a trivial consequence of the above definition. 

Theorem: Let rj and 6 be small quantities, and assume that A satisfies the thermody- 
namic normality (11. 7p with ( = 27]"^ 6. Then for any initial state tf{0) € 'He,ae, one can 

This terminology was suggested to the author by Joel Lebowitz and Sheldon Goldstein. The 
adjective "thermodynamic" stresses that we are dealing with a macroscopic observable in a macroscopic 
system where thermodynamic description is meaningful. We remark, however, that the thermodynamic 
normality is neither a sufficient nor a necessary condition for the system to exhibit thermodynamic 
behaviors. Also note that, in [7!, the term "normal" is used to indicate a more detailed property. 
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find a (sufficiently long) time T, and a "good" subset Q C [0,T] witlo \Q\/T > 1 — 5, 

{^{t),{A-{AU?v(t))<{Ar^f (1.9) 

holds for any t E Q. 

Roughly speaking, (11.91) says that for any t G ^, we have 

A = {A)^, + 0{A7]) in^{t). (1.10) 

Thus if one measures A at t, the outcome must be close to {A)^^ with probability very 
close to 1. Moreover the set Q of such "good" t occupies most of the time interval [0, T]. 
Therefore, the result of quantum mechanical measurement of A after a long enough 
time is almost certainly almost identical to the microcanonical prediction {A)^c- This 
justifies the use of the microcanonical averagj^. 

Note that the "good" set Q inevitably contains some intermittent vacancies. A 
perfect settlement to equilibrium is not possible because quantum dynamics is quasi- 
periodic. The present consideration does not provide any information on how large T 
should be to see the approach to equilibriuno. 

Let us be more precise about the above claim about the measurement. Let Pa be the 
projection operator onto the eigenspace corresponding to the eigenvalue a of A. Then 
Vaif) '■= {v{'t)y Pa, ^{t)) is the probability that the outcome of the measurement of A is 
equal to a. The inequality (II. 9p implies 

J2i^-{AUcfpa{t)<{Avf, (1.11) 
a 

and hence for any s > 0, one ha: 



Probt 



a-{A)^,\>As\:= ^"(^^-(l)' ^^-^2^ 



(\a-{A)n,c\>As) 



By taking s so that ?7 ^ s ^ 1, the right-hand side becomes small, which means that 
the (quantum mechanical) probability of observing a value far from (A)jnc is small. 



Proof of Theorem: From (II. 3p . we have 

a,a' 

\G\ denotes the Lebesgue measure (the total length) of Q. 

The present formulation automatically covers the canonical average if one regards the system as 
consisting of a subsystem and a heat bath. When the subsystem itself is macroscopic, a macroscopic 
quantity A of the subsystem exhibits negligible fluctuation. 

It may be possible that a completely isolated quantum system sometimes requires quite a long 
time to relax to the equilibrium. 

The argument used here is that for the Chebyshev's inequality. 
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Since the energy eigenvalues are nondegenerate, the long time average of f ll.lSp becomes 
hm i C dt{ip{t), {A - (i)^e)VW> = E (^"' - (^)mc)'V^«> < 2{Ar,f5, 

Ttoo J Jq 

(1.14) 

where we noted that |cq,P = 1, and used the thermodynamic normality (11.71) with 
C = 2rf 5. Since we have the bound (11.141) for T oo, there exists T > such that 

I j\t{v{t),{A- {A)„,,)\{t)) < {Ar^f5, (1.15) 

Define the "bad set", on which A deviates considerably from {A)^^^, by 

B:=[te [0, T] I (^(t), (i - (i)^e)' ^(t)> > [Ar^f]. (1.16) 
Note that Q = [0,T]\B, and = T — \Q\. Then we have 



dt{v{t),{A-{AU^r^{t))>- j^Jt{A^f = {A^r-^ = {A^f(l ^ 

(1.17) 

which, with (I1.15p . implies the desired bound on \Q\/T. I 

1.3 A result which holds for most initial states 

If one is satisfied with a result which holds for most initial state (y9(0) from 'He,aei 
it is enough to make much weaker assumptions than the thermodynamic normality. 
Although this is beside the main point of the present note, let us briefly describe some 
facts for completeness. 

By examining the above proof, one easily finds that the inequality in fll.l4p is the 
key for the subsequent estimates. Thus we readily find 

Proposition 1: We do not assume the thermodynamic normality. If the initial state 
V5(0) = Y.a ^ T-iE,AE is such that 



E |c,|2 (i - (i)^c)' i'o) < 2{A7]f 6 (1.18) 



then the conclusion of the main theorem holds. 

We can also show that the condition fll.l8p is indeed satisfied for a typical state in 
'He,ae by assuming that the observable A exhibits small fluctuation in the microcanon- 
ical ensemble as in (II. 8p . 

Proposition 2: We do not assume the thermodjTiamic normality. Suppose that the 
inequality (II. 8p is valid with C, = 2rf6v, where i/ is a small constant. If one sam- 
ples a normalized state V9(0) uniforml}@ from the space 'He,ae, then Lp{0) satisfies the 
condition (I1.18P with the probability not less than 1 — v. 

We do not claim that the uniform samphng is physicaUy reahstic. The present proposition simply 
tells us that there are many (^(0) which satisfies the condition (|1.18p . 
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Proof: Write I := {a \ ^ [E, E + AE]}, and denote by (■ ■ ■ )c the uniform average over 
(Ca)ag7 with G C and J2aei 1*^"!^ ~ Note that the symmetry imphes {c^Ca')c = 
Sa,a'/\I\, where |/| denotes the number of elements in /. 

Write for simphcity fa '■= {'ipa, (^)mc)^ i^a)- Then the assumption of the propo- 
sition reads 

(EJc„|V.)c = ^E/" = ((^- <2{A^r6u (1.19) 

By denoting x[true] = 1, x [false] = 0, we see 

{EaK\'fa)c>2{A7^r6{x[EM'fa>2{Ar^rS]), (1.20) 
which with fll.l9p implies 

Prob[E„ Kl'fa >2{Ar^)'S] := (x[E« Kl'fa >2{Ar^yS] < u, (1.21) 
which is the desired bound. I 



1.4 Relation to the result of Goldstein, Lebowitz, Mastrodonato, 
Tumulka, and Zanghi 

Let us briefly discuss the results of Goldstein, Lebowitz, Mastrodonato, Tumulka, and 
Zanghi [7], and how our main observation in section [Tl2] is related to it (indeed in a 
straightforward manner) . 

Goldstein, Lebowitz, Mastrodonato, Tumulka, and Zanghi start by assuming that 
there is a subspace Tieq C 'He,ae in which various macroscopic quantities take their 
equilibrium values. They assume that the dimensions of the spaces satisfy 

dim Tien 

^ \. (1.22) 



dim'HE,AE 

This assumption is consistent with the picture that overwhelming majority of the states 
in 'He,ae can be regarded as representing the equilibrium (recall Figured]). Let Peq be 
the orthogonal projection onto the equilibrium Hilbert space "Heq- 
For a given equilibrium Hilbert space "Heq, consider the condition 

{ipa, Peq ^a) ^ 1 for all a such that E^ e[E,E + AE], (1.23) 



which, in our terminology^^ should be called "thermodynamic normality." 

Then the results of Goldstein, Lebowitz, Mastrodonato, Tumulka, and Zanghi [7] 
can be summarized (in our own interpretation) as follows. 



As we remarked before in the footnote [131 the authors of [7] use the term "normal" in a different 
manner. 
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(A) If the equilibrium subspace Heq satisfies the thermodynamic normahty fll.23p . then 
any initial state ip{0) G 'He,ae will eventually evolve into equilibrium in the sense that 

(^(t),Peq^(t)>^l (1.24) 

holds for most t in a sufficiently long time interval. 

(B) Fix the Hamiltonian H , and choose the equilibrium space T-ieq in a random mannei@. 
Then with a probability very close to 1, the equilibrium subspace satisfies fll.23p . 
Thus the thermodynamic normality is a typical property. 

Note that, as the authors of [7] themselves warn, the random choice of Heq in (B) 
should not be taken too literally. This formulation does not mean that either the space 
?^eq or the Hamiltonian H is chosen according to certain stochastic rules in reality. 
The typicality (B) should be understood as a strong indication that the thermodynamic 
normality fll.23p is a common property which is shared by many pairs of ?^eq and H. 

The nontrivial part of this work is (B), and the derivation of (A) is rather straight- 
forwarci^. Nevertheless, we were very much impressed by the claim (A). It can be read 
as a concrete criterion that "Heq and H should satisfy to guarantee the approach to 
equilibrium for any initial state. Then a natural challenge is to take a "constructive" 
approach, i.e., to choose concrete (and hopefully nontrivial) models and show that the 
criterion is actually satisfied (or show that it is not satisfied). 

To treat the expectation value of the projection Peq, however, is not an easy task 
in general. Since the equilibrium space "Heq is determined by possible macroscopic 
observables, it might be easier to work directly with observables than the projection. If 
one pursues this line of thought, and reinterprets the above (A) in terms of a macroscopic 
operator, then one gets our Theorem 1 in an almost straightforward manner. 

In the next section, we present the first step of our "constructive" approach by 
proving the thermodynamic normality of some observables in three concrete models. In 
the example of thermal contact treated in section 12. 3[ we can treat the projection Pgq 
explicitly and show the property fll.23p . In the free fermion model treated in section 
on the other hand, we still do not see how to treat Peq while we are able to show the 
thermodynamic normality for certain observables. 

2 Examples 

We shall discuss three classes of examples where we can establish the thermodynamic 
normality of certain observables explicitly. Thus they provide rigorous examples where 
the approach to thermal equilibrium from any V5(0) G 'He,ae is proved rigorously. 

In the original work [7], Hcq is fixed, and H is chosen in a random manner. But mathematically 
speaking, these two formulations are the same. 

Of course, to find the correct claim to prove is never straightforward. 
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2.1 Independent spins under random magnetic field 

Let us start with a trivial example of independent spins under quenched random mag- 
netic field. In this model independent precession of each spin causes the "approach to 
equilibrium" for certain observables. Although everything is trivial, it may be a good 
idea to look at the simple (but genuinely quantum mechanical) mechanism that realizes 
the relaxation-like behavior. Interestingly, the same model also offers a counterexample 
to the thermodynamic normality. 

Consider a system of N independent spins with S = 1/2. We assume ^ 1. For 
u = x,y,z and j = 1,2, . . . , N, we denote by Sj"'' the spin operator in the t/-direction 
of the j-th spin. We consider the Hamiltonian 

N 

H = 'Y.^^Sf^ (2.1) 

i=i 

where the quenched random magnetic field hj (j = 1, . . . , A^) is a random quantity 
which are independent with each other and identically distributed according to a certain 
continuous probability distribution. Since the details of the distribution is irrelevant 
here, one may assume that each hj is drawn uniformly from the interval [—h, h]. 

For each j = 1, . . . , A^ and a = ±1, we denote by ipj the standard basis states of the 
j-th spin which satisfy 

^fVj = |^J. (2.2) 

Let cr := {aj)j=i^...^N be the multi- index with aj = ±1, and define the corresponding A^ 
spin states by 

^.:= (g) (2-3) 

j=l,...,N 

It is obvious that the state \E'o. is an eigenstate of H with the eigenvalue 

N 

E^:=-Y,hjaj. (2.4) 

i=i 

Since each hj is random, the energy eigenvalues (12. 4p are nondegenerate with probability 
one. 

Let us take the total spin in the x-direction 5'^^^'' = Sj^^ as a macroscopic 

observable. Then one readily finds that 

(VI/., si:i = 0, (VP., {si:ir = ^ (2.5) 

^ (x) 

for any configuration cr. Thus by choosing the typical magnitude of Sl^t as A^/2, we find 
that the thermodynamic normality (II. 7p is satisfied with (S'^t )mc = and ( = <^ 1. 
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Note that the range of the energy is arbitrary here. One can even start from the 
state 

"^(0)= (8) ' ^ (2-6) 

j=l,...,N 

which satisfies "^(0) = (A^/2)$(0), and conclude that | (<l>(t), 5^ ^(t)) | < for 
sufficiently large and typical t. As we have noted in the beginning, this "relaxation" is 
nothing but a trivial consequence of independent precession of each spins. 

It is interesting to see what happens if we take the total spin in the z-direction 
^ul = our macroscopic observable. Since the energy eigenvalue (12. 4p is 

the sum of (the sign factor times) the continuously distributed random quantities, it 
happens in general that two energy eigenvalues Ea- and E„i which are extremely close 
to each other have radically different configurations cr and cr' . As a consequence [TO] , 
which is the eigenvalue S^^l, shows a rather erratic behavior when viewed as a 

^ (z) 

function of the energy eigenvalue Efj. This implies that the macroscopic observable S^^l 
is not thermodynamically normal in any energy intervals. Indeed, the expectation value 
of S^^l is independent of time, and can never approach to equilibrium. 

This rather trivial example illustrates how the thermodjTiamic normality (as well 
as the tendency of approach to equilibrium) can be lost in a system with a quenched 
disorder. In [10], an interacting spin model with the Hamiltonian 

N N 

H = jY,SrS,+i + Y,h,a„ (2.7) 

i=i i=i 

where J is a constant and hj is uniformly distributed in [—h, h], was studied numerically. 
A systematic analysis suggests that the "localization" observed above for J = persists 
in a model with sufficiently small \J\. In this case, the thermodynamic normality is not 
satisfied for a general observables, and the model lacks the ability to relax to equilibrium 
by itself. For large enough | J|, the system enters the "delocalized" phase where the ther- 
modynamic normality may hold. See [10] and references therein for further discussions 
about the localization in many body quantum systems and its relation to the problem 
of the approach to equilibrium. 

2.2 Free fermions on a discrete chain 

The second example is the system of free fermions on a discrete chain of L sites 
with the standard hopping Hamiltonian (12. 8 p (with an extra phase factor 6). We show 
explicitly that the particle number A^^ and the energy He on a part of the lattice are 
normal when N/L and i/L are 0(1) and the system size L is large. With a further con- 
sideration on the degeneracy, we prove that both Ni and He approach their equilibrium 
values in the sense of the main theorem!^, provided that L is odd and the phase factor 

We thus allow any initial state from He.ae- Although it is clear that TLe.ae contains various 
nonequilibrium states, it is not easy to precisely state which states are in 'He,ae- 
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6 does not take one of a finite number of exceptional values. 

Note that this conclusion is consistent with the well-known fact that an ideahzed gas 
never exhibits a true thermahzation towards equihbrium. We are here only looking at 
rather crude observables like the partial particle number and the partial energy. These 
quantities can (and do) approach its equilibrium values with only ideal gas dynamics. 
If we look at more sophisticated observables (such as properly coarse grained velocity 
distribution), then they must fail to be thermodynamically normal. 

It is well known that a classical ideal g classical system of a macroscopic 

number of non-interacting particles) also exhibits a diffusion-like behavior. Suppose that 
all the particles are initially located in one end of a container. As each particle moves in 
the container with its own velocity, the positional distribution of the particles gradually 
spreads out, and in the long run one finds that particles are distributed almost uniformly 
in the container. The relaxation phenomena we see in this section is quite similar to this 
classical phenomenon. It must be noted, however, that in the classical model, one needs 
to exclude exceptional initial states (like the one in which all the particles have the same 
velocity) to have the relaxation behavior. There seems to be an essential difference with 
the following quantum version where we are allowed to start from any initial state. 

2.2.1 Model and energy eigenstates 

We consider a system of N spinless fermions on a finite periodic chain with L sites. We 
assume both L and N are large, and the density N/L is of 0(1). 

Let us denote a site on the chain as a; G {1, 2, . . . , L}. The Hamiltonian is 

1 ^ 

^ = 9 Ei^^' < ^-+1 + <+i (2-8) 

where and c* are the annihilation and the creation operators, respectively, of the 
fermion at site x G {1, 2, . . . , L}. They satisfy the standard canonical anticommutation 
relations 

{<, Cy} = 6^,y, {<, C*y} = 0, {C^, Cy} = (2.9) 

for any x,y ^ {1,2, where anticommutator is {A,B} := AB — BA. We use 

the periodic boundary condition c^+i = Ci. (12. 8p is the standard tight-binding hop- 
ping Hamiltonian except for the phase factor 9 G (0,7r/L), which is (rather artificially) 
introduced to reduce the symmetry of the system and avoid degeneracy. 

For n G {1,2, define the creation operator for the single-particle energy 

eigenstate as 

1 ^ 

* x=l 

which again satisfy the canonical anticommutation relations 

{a*, an'} = Sn,n', {<, </} = 0, {a„, an'} = 0, (2.11) 



I — — X 



c 



(2.10) 
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for any n,n' G {1, 2, ... , L}. The fact that a* creates a single-particle energy eigenstate 
is encoded in the commutation relation 

[^,<]=e„< (2.12) 

with the corresponding single-particle energy eigenvalue 

en := coi^y-j- + Oy (2.13) 

Fix the particle number N . Let <l>vac be the vacuum state with no particles on 
the chain, which satisfies $vac = for any x G {1,2,.. For each subset F C 
{1, 2, . . . , L} with |F| = A^, one has an N particle energy eigenstatj^ 

'fr= (n<)'^'vac, (2.14) 

ner 

with the energy eigenvalue 

Er = Y,^n. (2.15) 

ner 

The energy eigenvalues fl2.15p may be degenerate. It can be shown that by choosing 
L to be a prime number, the energy eigenvalues are nondegenerate for most choice of 
^ G (0, tt/L) except for a finite number of exceptional points. In this case we can apply 
our main theorem as it is (provided that the thermodynamic normality is guaranteed). 

It is of course somewhat absurd to limit ourselves to chains whose length is exactly 
equal to a prime number. In fact, in order to prove the desired conclusion about the 
approach to equilibrium, it is enough to assume that L is odd and choose suitable 6 
except for a finite number of exceptional points. See section 12.2.31 

Although we agree that these fine-tuning concerning the degeneracy may not sound 
quite physical, we do not regard this as a serious problem. Technically speaking, it is 
obvious that the degeneracy can always be lifted by using various "dirty" tricks, e.g., by 
adding extremely small random on-site potential energy. It is our optimistic hope that 
the degeneracy (in a clean system) is an accidental phenomenon in a non-interacting 
systems, and may be naturally avoided in fully interacting systems unless there are 
special symmetries. 

2.2.2 Thermodynamic normality 

Let us discuss the thermodynamic normality (11.71) for the following two observables. 
Fix i < L, and write v := i/L. We assume v = 0(1). Let 

e ^ e 

Nr.= Y,c:c,, Hr.= -Y.{e''c:c,+, + e-^'c:^,c,} (2.16) 

x=l x=l 

The product of a* is ordered according to the order of n. 
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be the partial number operator and the partial energy operator, respectively, in the 
interval {1, . . . ,£}. By explicit (and slightly nontrivial) calculations, whose details can 
be found in Appendix IA.21 we see that 

{^r,Ni^r) =vN, {^r, {Ne? ^r) = {v + {v - v"^) N + 0{e) (2.17) 

for any F. Thus one has (A^^)mc = v N for any E and AE, and 

($r, {Ne - (iV,)^c)' $r> = {v - v^) N + 0(£) = 0{L). (2.18) 

This means that A'^ is thermodynamically normal with respect to H with the choice 
Ne = N = 0{L) and C = 0{l/L) < 1. 
Similarly we can show that 

{^r.Hi^r) =vEr (2.19) 
($r, {Hi? $r> = {v Erf + {v - v') J^i^^? + 0{i) (2.20) 

for any F. Then in the microcanonical ensemble in the range [E, E + AE], one has 

{He)^, = vE + 0{AE), (2.21) 

and 

($r, {He - {He)mcf $r> = [v - v') J](e„)' + v [Er - Ef + O(^) 

ner 

< {v- v"^) Er + 0{i) = 0{L), (2.22) 

where we noted that0 |e„| < 1. Thus He is thermodynamically normal with He = E = 
0{L) and C = 0{1/L) < 1. 

2.2.3 System with general odd L 

Let us consider a system whose size L is a (large) odd integer, but not necessarily a 
prime number. In this case the energy eigenvalues Er may be degenerate. Then in the 
long time average f ll.l4p . off-diagonal terms ($r, (^ — (^)mc)^ '^'r') with F 7^ F' but 
Er = Er' do not vanish. 

However Lemma 2 in App endix I A . 1 1 implies that, for G (0, ir/N) except for a finite 
number of exceptional points, the degeneracy Er = Er' for F 7^ F' is only possible when 
F and F' differ in more than two elements (i.e., |F\F'| = |F'\F| > 2). We shall show in 
Appendix IA.2I that for such F and F' one has 

(<l>r,iV,$r'> = 0, (<l>r,(iV,)2$r,) = (2.23) 

This seemingly innocent manipulation relies essentially on the fact that the kinetic energy in the 
tight-binding model is bounded from above. We expect that in the continuum model, where e„ is not 
bounded, the quantity corresponding to He is not normal. This, we expect, is an artifact of an ideal 
gas. 
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and 

(<l>r,i^,$r'> = 0, ($r,(i^^)'$r'> = 0. (2.24) 

Therefore we see that the unwanted off-diagonal elements vanish, and the equality in 
fll.l4p is valid as it is. We have thus shown that both Ni and Hi exhibit the approach 
to equilibrium in the sense of our main theorem. 

2.3 Two identical system in a thermal contact 

The third example is a simple model of two identical systems in thermal contact. We 
show that the operator H'^^\ the Hamiltonian of one of the two systems, is thermody- 
namically normal when the "volume" of the system is large. Then the main theorem 
guarantees that H^^^ approaches its equilibrium value, even if one starts from a highly 
nonequilibrium initial state. Thus the well known fact that "two bodies in contact with 
each other will transfer heat until their temperatures are the same" has been derived 
only from quantum time evolution (although in a highly artificial model). 

We note that the present model is constructed so that our analysis can be carried out 
easily, and is hence extremely artificial. We nevertheless hope that it somehow mimics 
basic structures of energy eigenvalues in more realistic systems. 

Let us first describe a single system. The energy eigenvalue is n eo with = 1, 2, . . ., 
where eo > is a constant. The energy level with ne^ is fi„ fold degenerate. To mimic 
a healthy macroscopic system, the degeneracy is assumed to behave as 

a„~exp[ys„], (2.25) 

where V is the "volume" of the system, which will be made large. The "entropy" s„ > 
is increasing in n, and satisfies the strict concavity 

(2.26) 

for any n. An energy eigenstate with energy ne^ is written as "^^(^nj) where n = 1, 2, . . ., 
and j = 1, . . . , VLn- We denote the corresponding Hamiltonian as H^^\ Thus H'^^'* "^^(nj) ~ 

We prepare an exact copy of this system and denote its Hamiltonian and eigenstates 
as H^^^ and ipf^i'V respectively. Thus H^'^^ = n' eo V^^^) ,,•,)• 

The total Hamiltonian of the system is given bjO H = H^^^ + H^^^ + e i/int, where 
f^int is the interaction that we shall define. 

Let m > 2 be an integer, and let n,n' G {1, 2, . . .} be such that n + n' = m. Then 
any state ^[^^j) ® "^AI^/ j/) is an eigenstate of the noninteracting Hamiltonian H^'^'' + ff*-^-* 
with the eigenvalue meo. The degeneracy of this eigenvalue is given by 

= ^n^n' (2.27) 

n,n'>l 
{n+n'=m) 

"^^ Rigorously speaking, this equation should be -ff = H'^^'> (g) 1 + 1 (g) H^^'> + eHint- 
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We denote by Tim this f^m-dimensional eigenspace. The interaction Hi^t is designed so 
that to leave each T-L^ invariant, and mixes up the basis states within it. 

Although it is likely that almost any "generic" interaction i^int will do, we shall give 
a simple example of ifjnt to be concrete (and constructive). Fix m, and renumber the 
basis states V'ln j) ® ^{n' j') (with the restriction n + ra' = m) of T-Lm in an arbitrary 
mannei@ and call them tpk with k = 1, . . . , fi^- Then we define 

fl/2 if = 1 

I U otherwise. 

We do the same thing for all m to define i^int completely. 

The interaction f l2.28p lifts the degeneracy in T-Lm completely, and the normalized 
eigenstate is given by 



/ 2 

i^m,i := Y.MtT-^'') (2-29) 



k=l 

where ^ = 1, . . . , Vtm- The corresponding eigenvalue of the total Hamiltonian H is 

Em,t = meQ + e cosi-^ -j. (2.30) 

Thus, by taking e <C eo, we see that the eigenvalues of H are nondegenerate. This 
peculiar "band structure" of the eigenvalues of H makes the analysis easy. Although 
the band structure is just an artifact of the construction, we hope that the basic structure 
of the eigenstates fl2.29p mimics that in more realistic system^. 

Fix an even m, and consider the microcanonical ensemble for the energy interval 
[m eo — e.meQ + e\. The ensemble consists precisely of the sates in T-Lrn- 

Let us focus on the observable H'^^\ the energy of one of the two systems. Clearly its 
microcanonical average is {H)inc = rneQ. We shall show that iy*-^-* is thermodynamically 
normal when V is large. 

Let ra, n' be such that n + n' = m. The number of the basis states of Urn of the form 
V'(,lj)®V^(n'j') is ^n^n' — exp [V(s„ + s„' )] . Now writc n = (m/2)+r and n' = (m/2)-r. 
From the concavity f l2.26p . it follows that the quantity s„ + s„/ = S(rn/2)+r + ■S(m/2)-r 
attains its maximum at r = and decreases monotonically as r increases or decreases 
from C@. We thus find 

Q^m/2)+r^{m/2)^r ^^-^^^ ^2 31) 

^m/2 ^m/2 

for any |r| > 1, where Km '■= 2sm/2 — S(m/2)-i — ■5{m/2)+i > 0. This means that when V is 
large the space Hm mostly consists of the states '?/'(m/2j) ®'^{m/2j')' ^^^""^ eigenstate 



For example, one can use the lexicographic ordering in n, j, and j' . 

One can avoid too explicit band structures by considering (yet artificial) model in the line of [5]. 
^® Proof: (|2.26p implies s„ — s„_i > s„+i— s„. By repeatedly using this, one finds Sp—Sp^i > Sg+i — s, 
for any p < q. This means Sp + Sq > Sp-i + Sg+i, which justifies the claim. 
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fl2.29p is an almost democratic linear combination of all the basis states, one readily 
finds that 

(V'^,,,(^«)2V^^,,) = (meo)Ml + 0(e-^'''")}, (2.32) 

which means 

- (^«)^c)'V'm,^> = (meo)'O(e-^'-"). (2.33) 

Therefore the partial energy H^^^ is thermodynamically normal with the choice H^^^ = 
meo and C = 0(e"^'^'") < 1. 

We conclude that H^^'^ approaches its equilibrium value m eo when we start from an 
arbitrary initial state from Tim- Note that Tim contains such states where one system 
has much higher energy than the other. 

Finally let us remark how the analysis of Goldstein, Lebowitz, Mastrodonato, Tu- 
mulka, and Zanghi |7] can be carried out in the present example. If we identify T-Lm above 
with 'He,ae, the equilibrium Hilbert space Tieq of [7] is the one spanned by the states 

^(m/2,j) ®^(m/2j') ^^^^ = !,•••, Slucc dim "H^ = and dim?^eq = (fim/2)^ 

the relations f l2.27p and fl2.3ip justify ( ll.22p about the ratio of the dimensions of the 
spaces. The desired thermodynamic normality fll.23p of "Heq is then apparent from the 
expression fl2.29p of the energy eigenstate. 

3 Discussions 

In the present note we have shown that the thermodynamic normality of a macroscopic 
observable guarantees that the measured value of the observable approaches the corre- 
sponding equilibrium value when one starts from any initial state which have energy 
concentrated around a certain macroscopic value. This is a straightforward reinterpre- 
tation of the result of Goldstein, Lebowitz, Mastrodonato, Tumulka, and Zanghi [7]. We 
also presented three models in which the thermodynamic normality can be established 
for certain observables, thus providing concrete examples in which the "approach to 
equilibrium" is rigorously justified. 

Here we followed von Neumann [8] and Goldstein, Lebowitz, Mastrodonato, Tu- 
mulka, and Zanghi [7] in concentrating (except in section II. 3p on results which hold for 
any initial state from 'He,ae- That we allow basically an arbitrary initial state is a great 
advantage. Once we have a concrete model satisfying the necessary criteria, then we can 
justify the use of equilibrium statistical mechanics without resorting to any probabilistic 
(or measure theoretic) argument^. This makes the whole scenario extremely clear. 

But, to be honest, we are still not sure how essential this point is. It also sounds quite 
reasonable that results like those in section [L3] or in |H|6], which hold for "most" initial 
states from 'He,aei are sufficient from a "practical" point of view since the "chance" of 
preparing an inhibited initial state is negligible. 

We still have to restrict ourselves to "overwhelming majority of time" in the long interval [0,r]. 
But this has no conceptually difficulty like "probability" that the initial pure state falls into a specific 
subspace of the Hilbert space. 
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It must be seriously questioned whether the thermodynamic normahty holds for 
realistic macroscopic observables in a realistic macroscopic system which exhibits the 
approach to equilibriunj^. After all the thermodynamic normality is a very strong 
condition, and one might criticize that we are just assuming what we want. 

For the moment all that we can hope to do is to study further examples and learn 
more about the nature of energy eigenstates in interacting macroscopic systems. This 
is, of course, a formidably difficult task from the present standard of theoretical or 
mathematical physics. Numerical results for the Bose-Hubbard model in [TT] may be 
interpreted as an indication that the thermodynamic normality fails in a non-random 
interacting quantum many-body system. 

It may finally turn out that in general the strict thermodynamic normality does not 
hold in the sense that there exit certain exceptional energy eigenstates which consider- 
ably differ from the equilibrium state. Then we will have to be satisfied with results, 
like those in section 11.31 or in [HE], which hold for "most" initial states. 

It is also interesting to look for sufficient conditions for the thermodynamic normality, 
and examine their validity. For example, suppose that the Hamiltonian and hence the 
energy eigenstates are translation invariant, and the macroscopic observable A is written 
as A = Ylf=i where each cij is a translational copy of a local observable oi. Then the 
assumptions that, for any a such that Ea G [E,E + AE], the local expectation value 
satisfies 

{ipa,aii!a) ^ a, (3.1) 
where a is a constant, and the truncated correlation 

Clj Clj' i)o) - {ipa, CLj tpa) {i'a, i'a) (3.2) 

decays fast enough in the distance between the supports of dj and dj/ are enough to 
guarantee the thermodynamic normality of A. 

When the system exhibits a phase coexistence, the thermodynamic normality cannot 
be satisfied as it is. To treat such a situation, it seems that we need to consider a 
projection operator explicitly. Suppose, for simplicity, that two phases + and — coexist 
at a macroscopic energy E. Then the microcanonical average can be decomposed as 

(i)^c^r/(i)ic + (l-^)(^)mc, (3.3) 

where (■ ■ ■ are expectations in each phase, and r] G (0, 1) is a constant. We then define 
-P4 as the projection operator onto the subspace where A takes a value close to {A)^^ 
or (v4)~j,. P| is thus a sum of two projection^. Then our notion of thermodynamic 



As we have remarked at the end of section 12.11 a class of systems with quenched randomness 
provide examples which lack the ability to relax to equilibrium by themselves. Then the thermodynamic 
normality cannot be satisfied. Here we are concerned (or worried) with those examples which do have 
the ability to relax but do not satisfy the thermodynamic normality. 

It may be nontrivial to define a suitable projection operator in a general situation where, for 
example, spatial phase coexistence (separated by a domain wall) takes place. 
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normality should be replaced by 

{ipa, Pa^o) — 1 for any a such that G [E,E + AE]. (3.4) 

But such a condition involving a projection operator brings us back to the formulation 
of Goldstein, Lebowitz, Mastrodonato, Tumulka, and Zanghi [7]. See section [Lll 

A Some details of the free fermion model 

We supply some technical details necessary to get the results of section 12. 2[ 
A.l Proof of nondegeneracy 

Throughout the present subsection we make the 9 dependence of the energies explicit 
and write en{0) = cos(27r(n/L) + 9), and Er{9) = J2ner^n{9). We also define and use 
the complex constant ( := exp[(27r/L)i]. 

The first lemma states that the energy eigenvalues have no degeneracy for most 9 
provided that the system size L is a prime number. 

Lemma 1: Let L > 2 be a prime number, and fix any particle number N < L. Then 
the following is true for 9 G (0,7r/L) except for a finite number of points. For any 
F, F' C {1, 2, . . . , L} such that |F| = |F'| = and F ^ F', one has Eri9) Et>{9). 

Proof: For each F C {1,2,...,L}, let ziV) := EnerC- We shall show below that 
z{r) 7^ for any F such that |F| < L, and z(F) ^ z{r') for F and F' such that 
|F| = |F'| < L/2 and F T'. 

Let us prove the Lemma when < L/2. The case N > L/2 then follows by noting 
that Y^^^i en{9) = 0, and hence Er = --E{i,2,...,L}\r- 

Consider all the subsets F C {1, 2, . . . , L} with |F| = A^. From what we stated above, 
z{r) corresponding to these subsets are all distinct and nonzero. Denote by Pn the set 
of these (^) distinct and nonzero points in the complex plane. By noting that Er{9) = 
3?[z(F) e*^], we see that the energy eigenvalues of the N particle system is obtained by 
first rotating Pjy around the origin by 9, and then projecting the resulting sets to the 
real axis. Since there are only finite number of points, the projected images may overlap 
but only for finite "unlucky" values of 9. This proves the desired nondegeneracy. 

It remains to show the claim about the properties of ^(F). Here we shall make use of 
the classical result by Gauss known as "the irreducibility of the cyclotomic polynomials 
of prime index" (see, for example. Chapter 12, Section 3 of [T2] or Chapter 13, Section 2 
of [13] )• It implies that the — 1 complex numbers C, C^? ■ • •) C^^^ are rationally inde- 
pendent. More precisely, if "^^=1 C" = with integers mi, ... , mAr_i, one inevitably 
has mi = 7712 = • ■ ■ = rriN-i = 0. 

To show the first claim, we choose and fix hq G {1, 2, . . . , L}\F. Such no exists 
whenever |F| < L. Then z(T) = E^gp C"""" is rewritten as Yln=i™n("' with 
(mi, m2, . . . , m^v-i) 7^ (0, 0, . . . , 0). From the rational independence, we see z(T) 7^ 0. 
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To show the second claim, we choose and fix ng G {1, 2, . . . , Iv}\(r U F'). Such a uq 
exists since |r| = |r'| < L/2. Again {^(r) - z{T')} = Ener C"~"" - E„er' C"~"° 
is rewritten as J2n=i C" with (mi, m2, . . . , rriN-i) 7^ (0, 0, . . . , 0). From the rational 
independence, we see z(T) — z(T') 7^ 0. ■ 

The second lemma states the property of the energy eigenvalues used in section 12. 2. 31 

Lemma 2: When L is an odd positive integer, the following is true for 6 G (0,7r/L) 
except for finite number of points. For any n, n', m, m' & {1,2, ... , N} with n < n' and 

m < m' , and {n,n') 7^ {m,m'), one has 

en{e) + en'{e)^em{e) + ern'{e). (A.l) 



Proof: Note that the statement of Lemma 2 is that of Lemma 1 restricted to = 2. 
Thus it suffices to show that the set P2 consists of L{L — l)/2 distinct nonzero points. 
This is easily proved by the following geometric consideration (we suggest the reader 
to draw simple figures while reading the following). Since L is odd, it is trivial that 
+ 7^ 0. To show the distinctness, fix n, n' , and examine the relation + = 
u + V where u,v,E C such that \u\ = \v\ = 1. Since the simultaneous equations 
\u\ = 1, IC" + C" ~ ""I = 1 for M has at most two solutions (corresponding to the two 
intersection points of two circles), we find that only u,v satisfying the above conditions 
are {u,v) = (CC') or {u,v) = (C',C)- ■ 

A.2 Evaluation of the expectation values 



We describe calculations needed to get the results in sections 12.2.21 and I2.2.3[ Although 
the required techniques are standard, some results require nontrivial estimates. To help 
the readers who do not have experiences in many-body quantum systems, we also explain 
standard estimates in some detail. For a more systematic exposition, see, for example. 



For i < L and r;,G{l,2,...,L}, define the "incomplete" creation operator 



* x=l 



27rn ' 
I —— X 



(A.2) 



which corresponds to a* defined in (12.101) . From the basic anticommutation relations 
fra . we find 

{a^,hl} = ^=v, (A.3) 
for any n. It is important to note that, for n 7^ n' , the anticommutator 



exp 



■ 271 (n'-n) ■ 
I X 



(A.4) 



x=l 
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is nonvanishing because the sum in x is incomplete. 

Let us describe how to calculate the quantities involving N^. 
From fl2.9p . we get the commutation relations 

Kcy.<] = ^y,zCl (A.5) 

for any x,y,z & {1, 2, ... , L}. This relation and the definition fl2.16p implies that 

[iV,, <] = &:, [Ne,b:]=bl (A.6) 

which will be useful. 

Take F C {1,2,...,L} with |F| = A^, and write it as F = {rii, 722, . . . , njv} so 
that rij < rij^i. By using the definition fl2.14p of $r, the relation iV^^vac = 0, the 
commutation relation (lA.6p . and the relation {a*, 6*,} = 0, we find 









■ ■ ■ <J "^-vac 


N 




= J2< 








N 












N 













ner n/Gr\{n} 

where a{r;n) = (—1)-'^^ is the fermion sign. To calculate the expectation value, we 
rewrite $r similarly as 

$r = a(F;n)<( J] (A.8) 

n'er\{n.} 

for each n G F. Then we have 



ner n'er\{n} n'er\{n} 



ner n'er\{n} n'er\{n} 

n <')^vac, ( n <')^vac> 

ner „'gr\{n} n'er\{n} 

t;Ar, (A.9) 
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where we used flA.Sp . We have also noted that a„ ^nn'er\{n} '^n'^^^ 



(n.'enW <' ) «n vac = 0. 



To evaluate the expectation value of (iV^)^, we proceed as in (]A.7p to get 



n'eT\{n} 



n,n'er 
(?i<?i') 



n"gr\{n,n'} 



(A.IO) 



where a{T;n,n') = ±1 is an appropriate fermion sign. Then we have 

($r,(iV,)2$r)=^(a:( J] <')<fvac,&:( J] ^'j'^vac) 

ner ri'er\{n} n'er\{n} 

+ 2J2{«'{ n <") "^vac, ( n "^j'^vac) 



n,n'er ra"gr\{ra,?i'} 
(n<n,') 



n"er\{n,n'} 



+ n <") "^vac, an' an fc*' ( H <")*vac>. 



n.n'er n"gr\{n,?i'} 
(n<n') 



n"er\{n,n'} 



(A.ll) 



By using the anticommutation relations ( ]A.3I) . ( ]A.4p repeatedly we get 



an' an &* &*/ = f ^ - \Wn,n'\'^ ~ V 6*, 0^' + Wn,n' a^' + a„/ 6* 6*, a„. 



(A.12) 



Since the non-constant parts acting on ^nn"er\{n n'} '^n") "^vac give 0, we conclude 



Wn.n' 



n,n'<^T 
(n<n') 



(n<n') 



(A. 13) 

Since {v N)"^ + {v — v"^) N is what one gets from an independent distribution of particles, 
the final sum is a nontrivial correction that we shall now control. 

Fix n and n' such that n < n' . Note that the oscillating factor exp[z {27r(n' — n) / L} x\ 
in (1A.4I) has the wave length L/{n' — n). When we sum this oscillating factor over an 
interval which is an integer multiple of the wave length, the oscillation perfectly cancels 
out and the result vanishes. This leads us to a simple upper bound 



J^exp 

x=l 



2n(n'-n) 
i X 



< 



L 



n ~ n 



which means 



\Wn.n'\ < 



'n' — n)"^ 



(A.14) 



(A.15) 
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For 77,, n' such that L/{n' — n) < £, we use this upper bound. For the remaining n, n' 
such that L/{n' — n) > £, we use the trivial bound | exp[z {27r(n' — n)/L} x] | < 1 which 
gives 



\Wn,n'\^< i^jj =V\ (A.16) 

By replacing the range n, n' G F by n, G {1, 2, . . . , L}, we bound the correction term 
as 



n,n'eT n,n'G{l,2,...,L} 

(™<n') {n<n') 



Wn.n' 



;{1,2,...,L} n,n'G{l,2,. 
-n>L/e) {L/e>n'-r 



n,n'G{l,2,...,L} ^ ^ n,n'G{l,2,...,L} 

(n'-n>L/£) {L/^>n'-n>0) 

<L dx — + L-v^ = 2i (A.17) 

Substituting this back to flA.13|) . we get the desired estimate 

(<l>r, {Nef $r> = {v + {v - v^) N + 0(£). (A.18) 

Finally we briefly discuss the derivation of ( 12.231) about the off-diagonal matrix 
elements. Let F and F' be such that |F| = |F'| = A^, F 7^ F', and Er = Er'. Then as in 
(lA.lip . we get 

($P,(iV,f $r>=$^(( n <')'^>vac,a^&:( n <')^^^c) 
Tier n'er'\{m} n'er\{n} 

n,n'er n"£r'\{m,m'} n"er\{n,n'} 
(n<n') 

(A.19) 

where, in the first sum, m is chosen for each n so that m ^ F\{n}. Likewise, in the 
second sum, m, m' are chosen for each n, n' so that m, m' ^ r\{n, n'}. Then we use the 
anticommutation relations (lA.Sp . (]A.4p to re-order 0^,6* and am' amh*^})*^, as in (1A.12P 

to get 

(<l>r',(iV,)2$r)=5^const.(( J] "^'j'^vac, ( J] <)^^^) 

ner n'er'\{m} n'gr\{n} 

+ 5^ const.(( W <„)<l>vac, ( n ^'j-^-vac). (A.20) 



n,n'er n"er'\{m,m'} n"Gr\{n,?i'} 

(n<n') 
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But the right-hand side is vanishing since we know (from Lemma 2 in Appendix lA.ip 
that 

T\{n} ^ r'\{m}, r\{n, n'} T'\{m, m'} (A.21) 
for any choice of n, n' G F and m, m' G V . 

The calculation of the quantities involving proceeds in essentially the same man- 
ner, except for some minor complications. Let us make a few remarks about the differ- 
ences. 

From the basic commutation relations (lA.Sp . we find 



\ e-^^ c\ if X = 1 



which implies 



[^^^c:] = <(i(e^^c:_i + e-^c:+i) ifx = 2,...,£ (A.22) 
\e'^ c\ ifx = £ + l, 



[^,, <] = 6„ 6: + e-(^+'=) -I={e^M^+i)c,Vi - e^^ct } =: (A.23) 



where we wrote k = 2TTn/L. Comparing with the corresponding relation (]A.6|) for 



the partial number operator, we see that there appear extra boundary terms. Further 
commutation gives 

[Hi, dl] = e„ 6* + (finite number of terms) =: /*. (A. 24) 

One then repeats the same calculations as above, carefully replacing 6* by li* or /*, and 
computing appropriate commutation relations. 

It is a pleasure to thank Shelly Goldstein and Joel Lebowitz for useful discussions 
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